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Abstract 

We investigate quantum depinning of the one-dimensional (ID) commensurate charge- 
density wave (CDW) in the presence of one impurity theoretically. Quantum tunneling 
rate below but close to the threshold field is calculated at absolute zero temperature by 
use of the phase Hamiltonian within the WKB approximation. We show that the impurity 
can induce localized fluctuation and enhance the quantum depinning. The electric field 
dependence of the tunneling rate in the presence of the impurity is different from that in 
its absence. 

KEYWORDS: charge-density wave (CDW), quantum depinning, commensurability, impu- 
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1 Introduction 

The tunneling of macroscopic object, the so called Macroscopic Quantum Tunneling (MQT), is 
one of the most dramatic quantum effects. There has been many theoretical studies of MQT. [pi 
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Most of them are dealing with regular systems. In many actual cases, however, some kind of 
randomness exists and is expected to play an important role. Here we study the effects of 
randomness for the case of charge-density wave (CDW) as a typical example. |2j MQT of the 
commensurate CDW below but close to the threshold field was investigated by Nakaya and 
Hida (NH). || They employed the WKB approximation in order to obtain the tunneling rate 
as a function of the electric field. However the effects of randomness were not considered there. 
In the framework of classical mechanics it is shown that the randomness can make lower the 
threshold field for the depinning of the commensurate CDW with one impurity. [|J In this paper 
we calculate the quantum tunneling rate of the commensurate CDW with one impurity. The 
ID phase Hamiltonian J?] is adopted, and the electric field dependence of the tunneling 
rate is obtained by use of Langer's method || at absolute zero. It is found that the impurity 
can induce localized fluctuation and enhance the tunneling rate near the threshold field. We 
have different electric field dependence of the tunneling rate in the presence of the impurity 
from that given by NH in its absence. M 



2 The Model and the Threshold Field 

The CDW is the modulation of the charge density, p(X, T) = p cos (2kpX + (p(X, T)), where 
Po is the amplitude of the modulation, <p is its phase and kp is the Fermi wavenumber. We 
investigate the ID commensurate CDW with one impurity located at Xj at absolute zero 
temperature. The action of the system is described by <fi as 



S= dT dX 



A 



dT 




-g(l- cos(M0)) + V cos(2fc F X + 4>)5{X - X { 

The first term in the integrand is the kinetic energy density with A = (hv-p) / '{AnvV) , where 
t>F is the Fermi velocity and v<f> is the phason velocity. The second is the elastic one with 
B = (Hvf) I (47r) and the third is the electric field one. Here F e i e = (eE ext )/ (27r) and e < is the 
charge of an electron and E ext is the electric field. The fourth represents the energy resulting 
from the commensurability, M = 7i/(kpa), and a is a lattice constant. The last is the impurity 
potential energy with the amplitude V. In this action there exist the characteristic length and 



time, £ = y (2B) /(M 2 g) and £ t = J (2A) /(M 2 g), which are the phase coherence length and 
time due to the commensurability, respectively. We scale this action by £ and £ i; i.e. x = X/£, 
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t = T/£f The Euclidean action is introduced as Se = iS; 



S F = C / dr / dx 




+jH (! ~ cos(M0)) - vcos(x(xi) + 4>)8{x - x { 

where r = —it is the imaginary time, C = M 2 g££ t = 2\/AB, e = Fj (gM 2 ), and v = V/ (g£M 2 ). 
The typical magnitude of C is given as 



C/h = ^Jm*/m h /27T ~ 2 ~ 2 x 10. (3) 

Here m* and wib are the effective mass of the CDW and the band mass, respectively. We define 
Xi = Xi/£ and x( x i) = {2ttz)/M, where z = Xi/a characterizes the location of the impurity 
relative to the crest of the undeformed CDW with minimum commensurability energy. The 
range of x can be set as —tt/M < x < tt/M, because the action except of the impurity term 
has the periodicity of 27r/M with respect to the phase variable. We consider the case of v > 
hereafter. 

We first consider the threshold field at absolute zero temperature in the framework of 
classical static mechanics. Q The first variational equation of the action, 

- 0" - e + — sin(M0) + v sin(x + 4>)8{x - a*) = 0, (4) 

determines the ground state configuration of the phase with the boundary conditions, 0(±oo) = 
4>o and 0'(±oo) = 0. Here 

1 . . 

^° ^ M arcsm ( eM ) 

is the ground state configuration in the absence of the impurity. 

Depending on the range of Xi there are two kinds of configurations. If x is i n the range 
of — 7r/M < x < 0, the solution is located in the positive side (Config.l) as shown in Fig. [l| 
If x is i n the range of < x < tt/M, the solution is located in the negative side (Config.2). 
We take electric field e > 0. Under this condition, Config.l has more tendency to depin than 
Config.2. Therefore, we consider only the case of —n/M < x < 0, because we are interested in 
the lowering of the threshold field. 

The ground state solution, <f) s (x), in the presence of the impurity which is determined by 
eq. (HI) is obtained by the connection at the impurity site of two non-uniform solutions in the 
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absence of the impurity, (f>i(x) (Fig. ^|). The non- uniform solution, <pi(x), is given by 



M 



e ( — arcsin(eM) - ) + — [ <Jl - (eM) 2 - cos(M0, 



(5) 



which is derived by integration of 



- + — sin(M0) = O, 



(6) 



with the boundary condition, 0(±oo) = O an d 0'(±oo) = 0. Integration of eq. (§) is obtained 
only by numerical treatment. When we carry out the integration practically, we need the 
information of the phase at the impurity site obtained analytically as follows. The ground 
state solution, <fi s (x), is given by 



<p s (x) = 4>i(\x - Xi\ + c), 



(7) 



where c represents the distance between the center of (pi and Xi. Substituting eq. (|7|) into eq. ([|), 
we obtain 

- 2#(c) + v sin( X + 0,(c)) = 0. (8) 

By use of eq. 



eq. 



is rewritten as 



M 



e — arcsin(eM) - z (c) + — [Jl - {eMf - cos(M0,(c)) 



M 2 



+ 2 wsin (x + M c )) = °- 



(9) 



This equation determines the phase, <f>i(c) = (p s (xi), at the impurity site of the ground state 
solution, instead of c, when e is given. Substituting this <f) s (xi) into eq. @, we have the 
derivative at Xi, and then accomplish the numerical integration with the boundary conditions 
at X{ . 

The threshold field is determined as the field where the ground state solution becomes 
unstable. The eigenvalue equation of the fluctuations, 8<f>(x), around the ground state solution, 
4> s {x), is expressed as 

d 2 

+ VJx) 



dx 2 



5(p = X8(j), 



where 



V s (x) = cos (M(f> s (x)) + v cos (x + 4> s {x)) 8(x — Xj), 



(10) 



11' 
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and A is the eigenvalue. The threshold field is determined by the condition that the lowest 
A = 0. In the absence of the impurity, V s (x) becomes 



A g = VI - (sM) z . (12) 

Then the threshold is = 1/M. At the barrier of the potential energy density in the 
absence of the impurity, 

U{<j>) = -ecj> + ^(1 - cos(M0)), (13) 



disappears. The threshold field, e c , in the presence of the impurity [[Ll[] is lower than that in 
its absence under the condition of —ir/M < \ < —tt/(2M) and v > 0. When v is fixed, e c /e^ 
goes to 1 as x tends to —n/{2M). We also note that the threshold field goes to a finite value, 
£oo, even if v — > oo. For M = 4 and x — — tt/M, ^oo/^t — 0.725, which is the lowest threshold 
of the present model. 

First we consider the reason of the lowering of the threshold field in an aspect of the 
fluctuation. Due to the existence of the impurity, the potential term of the eigenvalue equation, 
V B (x), is non- uniform in space over the length scale of £ (Fig. f|(a)). We note that in the 
present dimensionless notation, £ is unity. We solve eq. ( |TUD numerically by transforming a 
differential equation into a difference equation (the Numeov method) and the shooting method 
on A„. [T2j |T3| As the result two bound solutions, &q(x) and $i(x), (Fig. Mb) and (c)), are 



obtained at most under the condition of —ir/M < x < —tc/(2M) and v > below but near the 
threshold field. Other solutions are non-bounded solutions whose eigenvalues form a continuous 
spectrum above A g . The eigenvalue of the bound state is smaller than the lowest eigenvalue, 
A g , of the continuous spectrum. Then applying the electric field, the lowest eigenvalue, Ao, 
of the bound state becomes zero before A g becomes zero. In this case the threshold field is 
lowered by the impurity than that in its absence. Here we note that in the absence of the 
impurity V s (x) is uniform in space and there exists the continuous spectrum only. In the region 
of — 7r/(2M) < Xi however, there is no bound state solution of eq. (fLOD, and the threshold is 
the same as that in the absence of the impurity. 

The bound state corresponds to the localized fluctuation over the length scale of £ around 
the impurity (Fig. |](b)). Owing to the localized fluctuation of the lowest eigenvalue, the 
CDW starts local sliding at the threshold field. The onset of a sliding in the absence of the 
impurity accompanies the spatially uniform and temporally continuous motion of the CDW. 
In the presence of the impurity, however, the sliding takes place around the impurity, i.e. the 
local sliding. Once the local sliding sets in, the CDW gains the electric field energy, which is 
transferred to the kinetic energy. Therefore, the sliding expands to the whole system and the 
whole CDW completely depins eventually, when dissipation is small enough. In the present 
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model dissipation is not taken into account. The issue of dissipation is discussed in the final 
section. 

Next we clarify the meaning of e c in an aspect of potential curves in the presence of the 
impurity. [|J In the absence of the impurity the phase is uniform in the ground state, and 
hence the potential energy curve can be expressed with respect to the uniform phase variable 
as eq. fllTf). With the impurity, however, the phase is not uniform. Apparently we should 
consider the potential curve in a functional space in this case. We can reduce the problem, 
however, to a one-particle problem by adopting the phase variable at the impurity site, (pi, as 
a kind of generalized collective coordinate. For each value of (pi, we can determine the lowest 
energy configuration; we solve eq. (|6|) under the boundary condition, (p(xi) = (pi an d 



(j>'(Xi 



1 , \ 1 



e ( — arcsin(£M) - <f>A + — ( y/1 - (eMf - cos(M0 i 



which is just the same as eq. (H) at the impurity site. Substituting this optimal configuration 
into the Euclidean Lagrangian, we obtain the potential energy, C/ imp (0j), as the function of (pi, 



f^imp ( 



-2V2 r d0 v /E o (0) - v cos( X + (14) 



where 



Eo(4>) = e arcsin(eM) -«/»)+ ± (\/l - (eMf - cos(M0)) . 



The value of (pi that minimizes U- irap {(pi) is (p s (xi). Then we introduce a variable, X = (pi — (p s (xi), 
and consider [7; mp as a function of X hereafter. The curves of Ui mp (X) is shown in Fig. ||. For 
each electric field there is 0i_ max , and for (pi > 0j_ max no such an optimal configuration exists. 
Therefore once the barrier disappears at e c , the sliding sets in because of the absence of any 
stable configuration, and then the CDW starts to depin. The value of 0i_ max is same as the 
maximum value of (pi(x), (p* (see Fig. §). 



3 Quantum Depinning in the Presence of an Impurity 

We derive the effective action based on the one-particle picture near the threshold field, e c , and 
calculate the tunneling rate in the presence of the impurity within the WKB approximation by 
Langer's method. || The tunneling rate is enlarged and its electric field dependence is changed 
by the impurity from that in its absence. 
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3.1 The Effective Action near the Threshold Field 



To investigate quantum tunneling problems theoretically by the semiclassical approximation, 
we must choose a suitable tunneling process, or in other words, the bounce solution. If we 
can derive the effective potential as a function of a single variable, the bounce solution is easily 
obtained by solving the equation of motion of a single degree of freedom. We call such a variable 
as the tunneling variable. In the present case the tunneling variable is the phase variable at 
the impurity site, 0» or X. The potential curve has been already obtained as a function of the 
tunneling variable in eq. (14). 

In order to derive the kinetic term, we follow the method used to derive the effective action 
in the Tomonaga-Luttinger model with a single barrier. [14] We introduce the partition function 
of the present model as 



T><f) exp 



(15) 



where Se[4>] is the Euclidean action which is given by eq. @. The phase, 0(r, x), is expressed 
by the sum of the static part, <fi s (x), which is the solution of eq. (|J), and the fluctuation part, 
S(j)(r,x); 

4>{t,x) = 4> s (x) + S(/)(t,x). (16) 

Substituting eq. (|1^) into the action and expanding it in powers of <50(r, x), we obtain the 
action up to the second order of £0(r, x) as Se — S + S 2 , where 

'd(j) s (xy ~ 



So = CI drdx 



Here 



So = C 



V(<f>) 







j drdx 


8(f)(T,x)^- 



v dx J 

d 2 d 2 



dr 2 dx 2 



+ V((f>.(x)) 



50(r, x) 



(17) 



-e6 + 



M 2 



cos(M0)) — v cos(x + 4>)8(x — Xi 



In order to integrate the phase, 50(t, x), under the condition, 50(r, Xi 
auxiliary fields, A(r), are introduced and the identity, 



IT 



■4>s(Xi) 



(18 
X(t) 



VX(t) J VA(t) exp 



i /drA(r) (6^T,Xi)-X(r)) 



is employed. Then the partition function is expressed as 



Z ~ exp(-S /h) J V8(j){r,x) J VX{t) J PA(r) 



x exp 



S 2 /h + iC/h / drA(r) (5<f){T, x^ - X{t)) 
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We diagonalize S2 by use of the eigenfunction; 



r, x) = °m,n exp {iuj m r) $ r 



a; 



where $ n (x) is the eigenfunction of eq. fliTf ) (see eqs. ( |TI| ) and (P~8|) , and note that V 
V s (x)), and c mj „ are the expansion coefficients. Then the partition function becomes 



(&(*)) 



Z ~ exp (-5 /n) / II dc m , n / VX{t) f PA(r) 



x exp 



1 C 

--r (| c m>n | 2 (u4 + A n ) - 22c m>n A(cj m )$ n (x i ) + i2A(u; m )X(-^ r , 



2 ft 



where 



A(u Tl 
X(u n 

Twice Gaussian integrations bring 



drA(r) exp(iu; m r). 
drX(r) exp(zco> m r) 



Z ~ exp / PI(t) / DA(r) 



x exp 



1 (7 E / A(^)A(-c m ) ! ^ 12 +z2AKOX( _^/ 



2 ft 



exp (— Sq/H) VX(r)exp 



ui m ) 

Here unimportant prefactors are neglected. Then the harmonic terms of the effective action is 
given by 

l^X(u m )X( 



5, 



~^E 



harm V- [gnfelP 

m 



(19) 

First we can now carry out the summation over n in eq. (0). Equation (|T9|) is expressed 
as a sum of the discrete eigenvalue part and the continuous one; i.e. 



$n(gj) I 2 



E 



$ (xi) | 2 L_ dfc | $(fc;a^ |2 



u& + A 



2tt 



a£ + A(fc) 



(20) 
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where L is the system size. Note that the first-excited state solution, $i(x), has odd parity 
and does not appear in the above equation. In eq. ( pOj) , A; is a continuous parameter which 
characterizes the continuous eigenvalues instead of n. The integration of the second term can 



be carried out analytically with the approximation, ~ yl/Lexp(ikxi). This is valid 

in the low temperature regime, which we are interested in, because the continuous states are 
higher energy states and give just a corrective contribution. Then 

l r a !*(*;«.) [■ i r- ik _ i 



2W-00 u^ + X(k) 2W-00 u^ + k 2 + A g 

1 1 



Next we consider the summation over uj m . In order to obtain the tunneling rate, we calculate 
the bounce solution and substitute it into the action. The bounce solution has a characteristic 
time scale, to^ 1 , which is shown in Fig. || Then the main contribution to the summation over 
oj m comes from terms with u m ~ u; B . Now we expand Sharm into powers of u m up to the second 
order; 

E |^ m)X 1 ( ~? ) - E(^(*& + K )X(u m )X(-u m ), (21) 



where 



*o(^) i 2 +-4 



(l $0(2*) I 2 + 



Ao 



and 



An 



$o(^) | 2 + 



This approximation is relevant under the condition of c^l/An 1- Here o>b is the characteristic 
energy of the motion of the phase variable at the impurity site, X(t), and Ao is that of the 
spatial fluctuation, $o(%), around the ground state. Then condition, c^l/An <C 1, means the 
spatial fluctuation around the ground state follows the motion of X(t) adiabatically. This 
condition is checked self-consistently later. 

Finally we obtain the total effective action as 



•SteLi — C / dr 



-M f— J + U imp (X) 
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where U] mp (X) is the potential energy given by eq. (|14]) with X = 0«— <j> a {xj). Here the harmonic 
term of the potential, KqX 2 in eq. (PH), is removed in order to avoid a double counting. We are 
interested in the tunneling near the threshold field, e c , where we can expand U imp into powers 
of X and obtain 

^imp(^0 — ^imp^ ^imp^ ; 

where 



"imp 



A {-e + jrsm(M(f) 



and 



Here 

Ec. is 



'imp 



6 



wsin(x + si ) 
4(-e+ jgsin(M0 si ))) 2 4cos(M0 



+ u cos(x + <Psi) 



+ 



+ wsin(x + 0si) 



(wsin(x + st ))' 3 ' vsin(x + s 
Xj), which is determined by eq. (§) for each field. Then the effective action near 

2 

2 " V dr 



cff 



C / dr 



— Mq [ — — I + CLimpX 2 — foimnJ^ 3 



Amp^ 



(22) 



3.2 Quantum Tunneling Rate near the Threshold Field 



We calculate the tunneling rate using the effective action near the threshold field. The tunneling 
rate within the WKB approximation is given by 

r imp oc exp [-S eS [X B )/h] 

where Xb(t) is the bounce solution of the equation of motion which is derived from S e g as 

d 2 x 



-Mn 



+ 1<x liavi X 3b[ mD X 



0. 



This equation of motion can be solved analytically and gives the solution as 

A B 



where 



-4, 



cosh 2 (ubt) 



"imp 



'imp 



1 "imp 

2Mq' 
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Substituting Xb(t) into «Seff[X] and carrying out the r-integration, we obtain the tunneling 
rate as 



Timp exp 



c 

Tl5 



-a- 



imp 



(23) 



The obtained tunneling rates are shown in Figs. |6] and [7[ Here we take the ratio of the effective 
mass of the CDW to the mass of the band electron as m* /m\> = 1 x 10 3 . This parameter 
determines a value of C (see eq. ([J])). The ^-dependence of the tunneling rate is shown in Fig. ^| 
for M = 4 and \ = —n/M, which gives the minimum value of the classical threshold field. On 
the other hand, the x-dependence of the tunneling rate is shown in Fig. [I] for M = 4 and v — 1. 
The condition, w^/Xq <C 1, that means the present effective action is relevant, is satisfied in 
all the region in Figs. ^| and |7]. We also note that the WKB approximation is relevant in the 
region of - lnr imp = S cS [X B ]/h < 1. 

The results shown in Figs. ||| and [7| can be expressed as 



r imp oc exp 



a(v,x)' 



£c(v,x). 



(24) 



Here a is a kind of generalized critical exponent and 7 is some factor. The v and ^-dependencies 
of a are shown in Figs. § and [|, respectively. When v exceeds 2 in the case of x — —tt/M 
(M = 4) and x exceeds — 77r/(8M) (M = 4) in the case of v = 1, either of the two conditions 
corresponding to the relevances of the present effective action, cu|/Ao 1, and of the WKB 
approximation gradually becomes invalid at any electric field. This means the present effective 
model is relevant in a limited region in the parameter space. Dependencies of 7 on v and x are 
shown in Fig. |l^ and Fig. |Tl|. In the relevant region, a depends on v strongly, but does not on 



X- 

In the absence of the impurity, an effective theory near the threshold field was given by 
Nakaya and Hida (NH) based on the scaling argument. The effective action of this theory is 
obtained by expanding U(<f)), eq. (0), at the inflection point, which yields This effective action 
is scaled as 



Si 



E-cff 



drdx 



~x 2 




~x 2 



+ 2 1 dx , 



where = (ai/&i)0 cf r, f = y/air, x = y/cfix, eff = - vr/(2M) + (2/M) Ml - e/e T )/2, 



1 — e/ex)/2, and b\ = (1/3!)M. . Therefore we obtain the electric field dependence of 
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the tunneling rate, r , without any knowledge about the bounce solution as 



) 



i 



T oc exp — B 1 



NH evaluated B numerically as 



B 



'0 - 



C6.2 x 10 
h M 2 



in the present notation. We compare the present result with NH's. In the limit of vanishing 
impurity potential, a in eq. (|4]) tends to unity which was given by NH in the absence of the 
impurity. On the other hand 7 in the present model, around 6 x 10 1 , is about three times larger 
than that given by NH (B ~ 2.0 x 10 1 for the present parameter.). 

We conclude that in the presence of the impurity a is larger than unity in the absence of 
the impurity. In the small limit of v , however, a becomes close to unity. Values of 7 given 
by the present calculation is about three times larger than that given by NH. Then present 
results in the small v limit shows the consistency with those of NH qualitatively, however, does 
not show exact coincidence quantitatively. The present effective model is applicable when the 
localized fluctuation can follow the motion of the phase at the impurity site, X(t), adiabatically. 
Then the origin of the ^-dependence of a seems to come from the existence of the localized 
fluctuation over the length scale of £ around the impurity, as in the case of the lowering of the 
classical threshold field. @] We also conclude that the impurity can enhance the tunneling rate 
near the threshold field, because a > 1 in that case. According to the scaling argument by 
NH, H the lower dimensional system has larger a: spatial dimension = 1, 2 and 3 have a =1, 
3/4 and 1/2, respectively. Then a > 1 means that the impurity brings the ID system to the 
lower- dimensional system effectively. 

4 Conclusion and Discussions 

We have investigated the quantum depinning of the commensurate CDW with one impurity 
at absolute zero temperature theoretically. The impurity causes the different electric field 
dependence of the tunneling rate from that in its absence and enhances the quantum depinning. 
In this issue, the fluctuation around the ground state plays an important role. It is described 
by the eigenvalue equation (0). In the absence of the impurity the eigenvalue equation has 
no bound state. On the other hand, in the presence of the impurity under the condition of 
— 7r/M < x < — 7r/(2M) the eigenvalue equation has a bound state, which means there is a 
spatially localized fluctuation around the impurity. This localized fluctuation which is induced 
by the impurity is the origin of the local sliding around the impurity and assists the depinning. 
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Finally we discuss future problems. The model which has been treated in the present paper 
offers some fundamental understandings of the effects of impurities on the depinning of the 
commensurate CDW. However, it is an idealized model. In actual cases, the following three 
points are important; dissipation, a finite density of impurities and three dimensionality. The 
dissipation is neglected in the present model. In the model once the local sliding occurs, the 
sliding expands to the whole system by the energy gain from the electric field. If a strong 
dissipation exists, however, the local sliding will be affected. Realistic systems contain not only 
one impurity but a finite density of impurities. However, the present result is expected to be 
applicable to the commensurate CDW with dilute but with macroscopic numbers of impurities 
where the inverse of the impurity density is larger than the phase coherence length, £. The 
depinning will be triggered by the local sliding at the optimal impurity site, x = —it/M. The 
effects of the high density of impurities is remained as a future problem. In one-dimensional 
systems, the impurity, a zero- dimensional object, has stronger effects than in higher dimensional 
systems. Generalization of the present model to three dimensional systems is left for future 
studies. 
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List of Figures 



Figure 1: A schematic picture of two kinds of the configurations, Config.l and 2, with respect 
to x- The potential in the absence of the impurity is U — —e<f) + (1 — cos0)/M 2 . 



Figure 2: The ground state configuration of the phase, <p s {x), in the presence of the impurity. 
The solid line is (f> s (x) for a choice of M — 4, v — 1, x — — tt/M and e = 0.1. 



Figure 3: The potential curves in the presence of an impurity for a choice of M = 4, x = — tt/M, 
v = 1 and e/e c = 0.90. The origin of the potential is set at 4>i = 4>s{xi) = and Ui mp ((j) s (xi)) = 0. 
The approximated U imv> in terms of X 2 and X 3 (X = 0« — s (xj)) is represented by the dotted 
line. For (pi > 0j_ max , there is no optimal configuration. 
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Figure 4: (a): A typical example of V s (x) for a choice of M = 4, \ — —n/M, v = 1 and 
e/e c = 0.95. (b) and (c): Normalized eigenfunctions of the discrete eigenvalue. There are only 
two bound solutions, $o(x) and $i(x). 

Figure 5: One bounce solution. Its typical time scale is u^ 1 . 

Figure 6: The tunneling rates in the presence of the impurity for M = 4 and x = —k/M. The 
electric field dependencies are different each other for various v values. The dotted line shows 
the result of Nakaya and Hida. 

Figure 7: The tunneling rates in the presence of the impurity for M = 4 and v — 1. The 
electric field dependencies are not so different each other for various x values. 

Figure 8: The f-dependence of a for M = 4 and x — — n/M. As v is getting smaller, a becomes 
smaller. When v — 0, a(0) — 1, which is given by Nakaya and Hida by the scaling argument. 
When v exceeds 2 in the case of x — — tt/M (M = 4), the present calculation gradually becomes 
unreliable. 

Figure 9: The x-dependence of a for M = 4 and v — 1. In this case, a does not depends on 
X strongly. When x exceeds — 7ir/(8M) (M = 4) in the case of v — 1, the present calculation 
gradually becomes unreliable. Note that the left most data point (open square) should be 
considered as semiquantitative one because of irrelevancy of the approximation we took (see 
text). 

Figure 10: The f-dependence of 7 for M = 4 and x — — vr/M. Values of 7 distribute around 
the value, 6 x 10 1 . 

Figure 11: The ^-dependence of 7 for M = 4 and v — 1. The value of 7 monotonically 
increases as x decreases. Note that the left most data point (open square) should be considered 
as semiquantitative one because of irrelevancy of the approximation we took (see text). 
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